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Abstract 

The connection between supersymmetric quantum mechanics and the Korteweg 
de Vries (KdV) equation is discussed, with particular emphasis on the KdV 
conservation laws. It is shown that supersymmetric quantum mechanics aids in 
the derivation of the conservation laws, and gives some insight into the Miura 
transformation that converts the KdV equation into the modified KdV equa- 
tion. The construction of the r-function by means of supersymmetric quantum 
mechanics is discussed. 



I INTRODUCTION 



It is well known that both the Korteweg-de Vries (KdV) equation [1] and supersym- 
metric quantum mechanics [2] have intimate connections to the inverse scattering 
problem [3-6]. In Ref. [4] it was shown that supersymmetric quantum mechanics can 
be used to construct reflectionless one-dimensional potentials with arbitrarily pre- 
scribed bound state energies. On the other hand, one can solve the KdV equation 
by means of the inverse scattering transform [1], wherein one associates a solution 
of KdV with a Schrodinger potential. Further links between supersymmetric quan- 
tum mechanics and the KdV equation have been found in connection with Backlund 
transformations [4] , which can be used either to add a soliton to a solution of KdV, or 
to add a bound state to a one-dimensional potential. It has been shown that super- 
symmetric quantum mechanics allows one to obtain an expression for the r- function 
that comes up in solving the KdV equation [5], and to obtain one- and two-soliton 
solutions for an infinite family of equations related to KdV [6]. In addition, it has 
been noted [7] that the change of variable connecting the modified KdV equation 
with the usual KdV equation comes up quite naturally in factorizing the Schrodinger 
equation. 

We have two parallel ways of adding a soliton to an already existing n-soliton 
solution. (1) We can perform a transformation within supersymmetric quantum me- 
chanics [3-6]. (2) We can construct a r-function for the n-soliton solution which, when 
acted upon by a suitable vertex operator, yields a r-function for the n + 1-soliton so- 
lution [8]. Our purpose in the present paper is to explore the relation between these 
two methods. 

In Sec. II we review the construction of one-dimensional potentials using super- 
symmetric quantum mechanics. In Sec. Ill we derive conservation laws for defor- 
mations of a potential that leave the bound state energies unchanged, while Sec. 
IV introduces a hierarchy of equations of the KdV type as a means of generating 
such deformations. In Sec. V, we introduce the r-function and the vertex operator, 
and present a method for constructing the n-soliton r-function using supersymmetric 
quantum mechanics. We conclude in Sec. VI. 



II SUPERSYMMETRIC QUANTUM MECHAN- 



In this section we consider the construction of reflectionless one-dimensional potentials 
by means of supersymmetric quantum mechanics. We begin by writing the one- 
dimensional Schrodinger Hamiltonian 
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(1) 



in factorized form as 



H+ = A* A, 



(2) 
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where A is given by 

A = -l + m. <3) 

In order that Eq. @ hold, f(x) must obey the Ricatti equation 

f + f = V + . (4) 
The form of Eq. (|2]) suggests the introduction of a 'partner' Hamiltonian if_ given 

H_ = ~ + V-(x)=AAi. (5) 
ax 1 

The potentials V± are given in terms of / by V± = f 2 ± /'. 

It turns out that there are certain interesting relationships between the spectra of 
H+ and H . For suppose that ip+ is an eigenfunction H + . We then have 

A*Aif)+ = E + ^ + . (6) 

Multiplying on the left by A gives 

Ad[Ai/> + ] = E+[Ai/> + ]. (7) 

From this it follows either that ip- = Aip + is an eigenfunction of if_ with eigenvalue 
or that ■?/>- = 0. Now the latter of these two possibilities implies that 

o = <v_IV-> = W+l^+> = (^+\A^+) = £+<^ + >, (8) 

so that ■?/>- = only when E + = 0. It follows that H + and if_ have the same 
spectra, apart from the single E = eigenvalue which is present in the spectrum of 
H + , but absent from the spectrum of A similar argument shows that if ip- is an 
eigenfunction of then ip + = A'tp- is an eigenfunction of H + . 

The factorization method can be used to construct reflectionless potentials pos- 
sessing an arbitrary spectrum of bound states. To see this, suppose we have a Hamil- 
tonian with potential having n eigenvalues E = E\, E2, ■ ■ ■ , E n . We may 
construct from a potential having the n eigenvalues Ek, plus one more 
eigenvalue E n+ i. To do this, first choose the zero of energy such that V {1) -> 
as x — > ±00, and define V_ by = — E n+ \. We factorize the corresponding 
Schrodinger equation as above: define f(x) by 

f 2 -f' = V. } (9) 

and construct the partner potential 

V + = f + f. (10) 

The potential V_ has the n eigenvalues E% — E n+ \, E 2 — E n+ i, . . . , E n — E n+ \, while 
V+ has n + 1 eigenvalues consisting of the n eigenvalues of V_, plus the additional 
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eigenvalue E = 0. Denning = V+ + E n+ i yields the desired potential, possessing 
the n + 1 eigenvalues E\, E 2 , . . . , E n+ \. In this way, one can start from the constant 
potential = and by iteration build up a potential possessing an arbitrary 
spectrum of eigenvalues. 

As an example of this procedure we construct a potential possessing a single bound 
state with energy E = — k 2 . We begin from = 0, so that V_ = k 2 , and / obeys 

f-f = K 2 . (11) 
This may be linearized by the substitution / = —w'/w, yielding 

/ = — k tanh k(x — x ), (12) 

so that V + is given by 



V + = K 2 



1 — 2sech 2 /t(x — xq) 



(13) 



and the desired potential, possessing a single bound state, is given by 

V {2) = -2/t 2 sech 2 /t(x - x ). (14) 

To see that is reflect ionless, consider a plane wave solution for the potential 
V-, ip-(x) = e lkx . The corresponding solution for the potential V+ is given by ip + = 
Aty-. From Eq. fljjD, it follows that ip + behaves asymptotically like (ik =p n)e lkx as 
x — > ±oo. Since no term proportional to e~ jfe:c arises, the potential is reflectionless. 
Defining the reflection and transmission coefficients R and T by 

s — ► — oo; 

(15) 

x — > +oo, 
we find that R = 0, while 

By repeating the above procedure, it is possible to construct potentials possessing 
arbitrarily many bound states. It is shown in the Appendix that such potentials are 
reflectionless, and that furthermore they have transmission coefficients given by 

N ih — k- 

r( ^n-^, (it) 

where the N bound states have energies Ei = —k 2 . 
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Ill CONSERVATION LAWS FOR ISOSPECTRAL 
DEFORMATIONS OF A POTENTIAL 



Consider again the one-dimensional Schrodinger equation, with a reflectionless po- 
tential V(x) = —u(x): 

?p xx + [k 2 + u(x)]ij = 0. (18) 

We wish to derive certain quantities which will be conserved under any deformation of 
the potential which leaves the spectrum of bound states unchanged. These conserved 
quantites may be obtained via an asymptotic expansion of the transmission coefficient 
for large values of k. We write the wavefunction as 



ijj(x) 



exp 



ikx + 



(x') dx' 



and expand (ft as 



Substitution of (|T5|) into 



4>{x) = J2 



fn{x) 



a (2^) n ' 
yields the recursion relation 



n-l 



,/nH 



~fn,x ^ t fkfn—ki 



(19) 



(20) 



(21) 



k=i 



for n > 1, while f\ = —u(x). The first few terms in the series are 

h = -U, f 2 =U x , /; " 2 1 - ' 



u xx - u , f 4 = (u xx + 2u ) x , and 



,2 

h = ~{uxx + 3u 2 ) xx + ul - 2m 3 . 



Taking the limit x — ► ±oo in Eq. (|19|) and comparing with Eq. flT5|) shows that the 
transmission coefficient given in Eq. QT7D may be written in two ways: 



00 l roo N - 

logT(fc) = J2 T^TSZ / fnip) dx = J2 



logfi/c - Ki) - \og(ik + Ki 



(22) 



Expanding the quantity on the far right in a power series in 1/k and equating coeffi- 
cients of l/k n gives 



f n (x) dx 



if n is even; 

-^E£r< if n is odd. 



(23) 



The integrals vanish for even n because the integrands in this case are total derivatives 
of a function having the same values at x — ±oo, as we see for n = 2 and n = 4 
in the list above. It follows from Eq. (|2~3D that the integrals of the functions /„ are 
constants of the motion for any spectrum-preserving deformation of the potential 
u(x): since the eigenvalues k 2 are constant for such a deformation, the integrals must 



4 



also be constant. These integrals are precisely the KdV Hamiltonians, up to an overall 
multiplicative factor. The choice of this multiplicative factor is simply a matter of 
convenience, so we choose to define the Hamiltonians as 

1 f°° 

#2ti+1 p] = — 22n+l 

/ f2n+3(x)dx. (24) 

A specific form for f n (x) may be written for the single-soliton solution by noting 
that for a soliton with x = 0, 

= Ah ikx ~ ( tk - Kt&nhKX ) e ikx 
\ ik + k J 

= exp(ikx + \og[ik — K,tanh.K,x] — \og[ik + k]) = exp[ikx + / 4>{x')dx'} 

J oo 

Expanding 

00 f (x) 

(f)(x) = = (d/dx)(\og[ik - /ttanh ra] - \og[ik + «]) (25) 



gives the result 



/n = _(M!i_tanh" ra . (26) 
n ax 



IV EVOLUTION EQUATIONS 

We can use the results of the previous two sections to derive a set of spectrum pre- 
serving evolution equations for the potential u(x). We know from the last section 
that any such evolution equation must conserve the transmission coefficient T(k), 
and, consequently, must also conserve each of the Hamiltonians H 2n +i- A natural 
choice for an evolution equation conserving a specific Hamiltonian H 2 k+i is a member 
of the Korteweg-de Vries hierarchy, which may be written in the form 

Ut = Tx^r- (27) 

We have introduced the variational derivative 5/Su, defined by 

r°° 5H 

SH[u] = / —8udx, (28) 

J-oo OU 

so that if H[u] = /f^ h(u, u x , . . .) dx, 

8H-(d_^_dh_ 
5u f- u \ dx) i)u s -- 
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where u^ is the k th derivative of u(x,t) with respect to x. Eq. ( p7|) automatically 
conserves the Hamiltonian H 2 k+i since 



(it J-oo 5m ' J-oo 2 V 5m 

The proof that Eq. (^7]) also conserves all of the Hamiltonians is somewhat more 
subtle. 

To show that Eq. (^) conserves all of the Hamiltonians, we first derive a new 
expression for the transmission coefficient T(k). We again employ a factorization 
method to write the Schrodinger equation in the form 

<h - [k 2 + u])ip = AAty = 0, (31) 



dx 2 

where A = — -^ + /, and 

/, - f ~ k 2 = u. (32) 
Since AA^ip = implies A*if) = 0, we find that 

/ = --^logV. (33) 
ax 

Since ip ~^ T(k) exp(ikx) as x — > oo, we conclude that the transmission coefficient 
may be written as 

/oo 
(life + /) dx. (34) 
-oo 

It follows that if the evolution equation for / corresponding to Eq. fl2"T|) has the form 
ft = dJ/dx, then T(k), and all of the Hamiltonians, will be conserved. It turns out 
that this is in fact the case. To see this, we make use of the following observations: 
First, we note that the variational derivatives of the Hamiltonians obey the recursion 
relation [1] 



F±i = - — + 2u x + Au— — ^ = -M — ; 



2k-l 



dx 5u 4 \dx 3 dx J 5u 4 5u 

Next, we observe that the linear operator M may be written as 



(35) 



( dx "0 dx ( dx "0 dx ^ ^ 

Finally, note that the variational derivatives 5H/5u and 5H/5f are related in a simple 
way: we have 

r°° 5H r°° 5H / d \ 

SH = / — 5udx= / —( — -2f)5fdx = 

J-oo ou J-oo ou V dx J 

r°° ( d \ 8H r 00 5H 

- L(-ai- 2 ^ dx -Ljf Sfdx ' (37) 
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so that 

S -(-!-*)£ 

We can now derive an evolution equation for f(x,t). Beginning from Eq. (ETH , we 
find 

d A , d 5H 2k +i 



dx Su 



I ( JL _ 2 f) — (-— - 2f) 6H2k ~ l - k 2 9 5H2k ~ 1 
4 V dx ) dx \ dx J Su dx Su 

9 nr \ d ( l5H 2k -i , k 2 5H 2k -3\ , f4 d 5H 2k -3 



dx J dx\ 4 of A Sf J dx Su 

We have made repeated use of Eq. to lower the index k, and Eq. to express 
6/ Sum terms of S/Sf. We may continue lowering k in this manner until the "remain- 
der" term (i.e. the term which is not multiplied by (d/dx — 2f) on the right hand 
side of Eq. (|39|)) has the form d / dx(SH_i/ Su). At this point, the remainder vanishes 
and the sequence terminates. Hence the evolution equation for / has the form 

d „„\/„ 1 d A, .v, +1 2l 5H 2 ( k -i)-i 



An argument due to Miura, Kruskal and Gardner [9] can now be used to show that 
Eq. p0| ) implies 

/-IGB-^fe). ,4!) 

The argument runs as follows: We know from Sec. Ill that /, which is the logarithmic 
derivative of the wavefunction, may be expanded in an asymptotic series of the form 

f = ik + y . (42) 

h ( 2 ^) n 



Consequently, the quantity 



1 9 l+lu2l SH 2 ( k -i)-i 



may also be written as a series in the form 



A' 



F= Yl b n {x,t){2ik) n . (44) 



Substituting these expansions into Eq. (|40|) and equating the coefficients of each power 
of k to zero, we find that the coefficient of k N+1 is simply b^, so that &at = 0. But 
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this in turn implies that b^-i = = &at-3 = . . . = 0. As a result, we conclude 
that F = 0, and Eq. (|T|) follows. 

The evolution equation fl4l|) is of the form 

/. = (45) 

From Eq. fl34|) it then follows that the transmission coefficient is a constant of the 
motion. Since the coefficients of the asymptotic expansion ( p2|) of logT(fc) are just 
constant multiples of the KdV Hamiltonians, it follows that each evolution equation 
of the form ( |27| ) conserves all of the Hamiltonians. Consequently, each such equation 
also preserves the eigenvalue spectrum of the potential u(x). 

The above formalism permits a simple derivation of the time dependence of the 
one soliton solution of KdV under the evolution equations (p7|). The fact that each of 
these evolution equations conserves all of the Hamiltonians implies that the various 
evolutions commute with one another [1]. As a result, the function u may be thought 
of as depending on an infinite number of time variables, with the dependence on the 
various times fixed by 

Ut ^ = Tx^nT- (46) 
We illustrate this by considering the one soliton solution of the KdV hierarchy. We 
make the ansatz u(x, ti,ts, . . .) = 2/t 2 sech 2 (/tx + X)o° a 2i+ihi+i), and subtitute into 
the evolution equations (fj6|) to determine the constants aj. We find that Eq. fl46|) 
with k = gives, after some manipulation, 

d °° d 00 

2/tai— sech 2 (/tx + c^m^+i) = 2/t 2 — sech 2 (KX + ^ (47) 
ox l=0 ox l=0 

so that «i = k. Using the recursion relation (|35D , it can be shown that 

~ k 2 - (48) 



dx 5u dx 5u 

From this all of the constants ctj can be determined from a\. The result is a.21+1 = 
k 21+1 , so that the full time dependence of the one-soliton solution is given by 

u(x, tx, t 3 , . . .) = 2/€ 2 sech 2 (/tx + Kti + k% + n 5 t 5 + ...). (49) 

An equivalent result has been given in [6]. 



V THE r-FUNCTION AND THE VERTEX OP- 
ERATOR 

In this section we discuss the construction of reflectionless potentials and multisoliton 
solutions of KdV using the r-function. A generic potential can be expressed in terms 
of the r-function by 

V{x) = -2— \ogr{x). (50) 



S 



The T-function can be constructed in either of two ways: First, we can use supersym- 
metric quantum mechanics to build up the r-function from a sort of "vacuum state" 
using the techniques discussed above. Alternatively, we can use the vertex operator 
[1,8] to construct the desired r-function. 

To build up the r-function using supersymmetric quantum mechanics, we consider 
a sequence of potentials Vi, V 2 , . . . V n , where V m has m bound states, the highest m — 1 
of which are shared with Kn_i. From the appendix, these potentials may be written 
in terms of a set of functions f m as 

Vm{ x ) = fm + fm ~ K m (51) 

or 

V m (x) = f m+ i — fm+l ~ K m+V (52) 

Using this, V n can be expressed as 

n d n 

V n =J2 ( V ™ ~ Vm-l) = 2— ^ f m , (53) 
m=l aX m=l 

where Vq = 0. Introducing functions w m defined by 

f m = -w' m /w m , (54) 

the potential V n becomes 

dx 

so that the r-function corresponding to V n is given by 



d 2 

V n =-2— loglJ^m, (55) 

m=l 



n 



II W m- (56) 



m=l 



We can use the operators A m = —dj dx+f m to construct the functions w m from simple 
linear combinations of exponentials. First observe that the sequence of Hamiltonians 
H m = —d 2 /dx 2 + V m may be written as 

H m — ^m-^-m — K m — A m+ iA m+ i — K m+ i- (57) 
Furthermore, the functions w m obey 

[^m^m — K m\ W m = ~K m W m - (58) 

Now suppose is chosen to satisfy 



d 2 

—w° m = [AtAl - K 2 ]w° m = -K 2 m w° m . (59) 



Multiplying on the left by A\ and using Eq. fl5"TD , we find that obeys 

[A 2 At - K 2 )A{w° m = -K 2 m A{w° m . (60) 
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Repeating this process we find that the function w m is given by 

n-l 

w m = n Aiw° m , 

k=l 

so that the r-function can be written as 

n /m— 1 



(61) 



n mm 



t „„o 



k l W m- 



(62) 



m=l v fc=l 

Since the functions are solutions of Eq. (0), they have the form 

< = a+e^ + a-e-^, (63) 

with the constants chosen so as to ensure that w m has no nodes. For the case of 
potentials symmetric about the origin, we choose a+ = 1, a~ = (— l) m+1 . The first 
few r-functions are 

r (x) = 1, (64) 
n{x) =a+e KlX + a~e- KlX , (65) 



and 



t 2 (x) 



r 3(£) 



+ a x a 2 (n 2 + Ki)e 



Ki)e" ± ~e'^~ + a{a 2 (—k 2 

RlX e K2X + a^a 2 (-K 2 + Kx)e~ Klx 



= a+a+a+(K 3 - k 2 )(k 3 - «x)(« 2 - Kl ) e ^+^+^ 

+ a+a+a 3 -(-/t 3 - Ka)(-«3 - «i)(«2 - Kl )e KlX+K2X ~ K3X + 

+ a+a^a^ (/t 3 + k 2 )(k3 - /«i)(-/«2 - ti)e 
+ 



KIX — K2X+K3X 



aja 2 a 3 (-« 3 + k 2 )(-k 3 - Ki)(-k 2 - /«i)e Ki:r ~ 



+ 



+ 



+ 



+ a x a^" a 3 ( K 3 - «2)(«3 + K i)( K 2 + «i)e 

+ ara+aj (-k 3 - k 2 )(-k 3 + Ki)(k 2 + Kl ) e -«is+~2*-K3s + 

+ a^a^ + « 2 )(k 3 + «i)(-«2 + Kl )e~ Kll - K ^ +K ^ + 

+ a^a 2 "a 3 "(-K 3 + « 2 )(-K3 + Ri)(-k 2 + K^e - " 1 ™ 3 ^ 3 *. 



(66) 



(67) 



The inductive generalization to the n-bound state r-function is evidently (cf. [10]) 



171=1 



Q t ^ rn ^m-£ 



m—1 

n 

Z=l 



5 rn. ^r. 



(6* 



An alternative method of adding a soliton to a r-function is provided by the 
vertex operator, which is a linear operator that, when applied to a r-function, adds 
a soliton. In Eqs. 



p7\) we have displayed the first few r-functions generated by 
super symmetric quantum mechanics. The given expressions are somewhat deceptive, 
however, since they over-count the number of degrees of freedom that one has when 
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adding a bound state to a potential. From the given expressions, it would seem that 
one is free to independently choose all three of the constants a+ , a~ , and K m . Contrary 
to this, one in fact has only two degrees of freedom: since the multiplication of a r- 
function by an overall constant leaves the potential unchanged, the true independent 
variables are the ratio a^/a^ and the constant K m . The former of these two is in fact 
a measure of the "position" of the soliton that one adds when adding a bound state. 
With these facts in mind, we re-write the above r-functions in an alternative form: 



Tl (x) = 1 + e 2Kl( - x ~ Xl) , (69) 

\K 2 + Kx. 



T 2 (x) = 1 + e 2K ^ x ~ x ^ -|_ e ^2{x-x 2 ) _|_ (^2l 1l\ e 2K 1 (x-x 1 )+2K 2 (x-x 2 ) ^ 



and, 



r J x \ = l _|_ e 2Ki(x-Xi) _|_ e 2K 2 (x-x 2 ) _|_ e 2K :i (x-x 3 ) _|_ f K 2 K l \ e 2n 1 (x-x 1 )+2K 2 (x-x 2 ) 



\K 2 + K 



_|_ ( ^3 ^\ 2 e 2 Kl ( a; -x 1 )+2 K 3( : r-a ; 3) _|_ ( ^ ^V^Cx-sa^Ks^-a*) 



+ / /t 2 - KA 2 / K 3 - Ki \ 2 / K 3 - ^2 \ 2 c 2 Kl (x-q :i )+2 K2 ( a: -x 2 )+2 K3 (x-X3)_ /y-Q 



In deriving these expressions from those given above, we have made use of the fact 
that multiplication of a r-function by an overall constant or by an exponential whose 
argument is linear in x leaves the potential unchanged. Using these expressions, we 
can construct a linear operator that converts an n-soliton r-function into an n + 1- 
soliton r-function. We denote this operator by A(k, x) to make its dependence on the 
relevant variables explicit, and construct it so as to satisfy 

A(K,x)r n = r n+1 . (72) 

Certain properties of A(k,x) can be deduced by considering its action on To = 1: we 
require that 

A(k 1 ,x 1 )t = A(K 1 ,x 1 )l = l + e 2 ^ x ~ Xl \ (73) 

This equation is automatically satisfied by A of the form A(k, x) — 1 + B(k, x), with 
B constructed such that 

£(«i,zi)l = e 2Kl(x - Xl) . (74) 

Further constraints are obtained by considering the action of A(k,x) on t\. Using 
(|74]), and applying A to n, we find 

A(«a, x 2 )t x = 1 + e 2Kl{x ~ Xl) + e 2K2{x ~ X2) + B{k 2 , x 2 )e 2Kl ^ Xl) . (75) 

B must therefore obey 

B(K2, X 2 )e 2Kl{X - Xl) = ( K2 ~ Kl ) 2 e 2 K1 (x~x l)+ 2 K2 (x-x 2 )^ (76) 

\K 2 + KiJ 
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Here we encounter the difficulty that such a constraint on B leads to inconsistencies 
when we apply it to r-functions of higher order. A way around this difficulty [11] is 
suggested by re-expressing the K-dependent coefficient in r 2 in the following form: 

K 2-Kx \ 2 2 (i og{ i_ Kl/K2 )_ lo g ( i +Kl/K2 )) -4Er=0 2WT(^) 2fc4 



K 2 + K\ 



e *Zjt=o 2i;+iU 2 / _ (77) 



We can generate precisely such a factor through the action of B by introducing a 
set of auxiliary variables £3, £5, £7, ... in the r- function according to the following 
prescription: replace the constant of integration x\ in t\ with the "time-dependent" 
expression Xi(k, £3, £5, . . .) = Xi — K i fc £2fc+i9 with Xi constant. These auxiliary 
variables can then be set to zero in the final expression for the potential. The k 
dependent constant can be generated by making use of the identity [based on Eq. (|77|)] 

~ 2 (7T 2 -£ + 5^|sfj + i^|s%+-) 2«i(a-a!i(/ti,ti ) ...)) = ( K<2 ~ Kl \ 2 ^^(x-x^^fy,...)) 



This is almost what we need. The correct form of B is evidently 

= e 2 (^fc°=o K2fc+1 * 2fc+1 ~ K5i )e -2 ^ fc=0 (afc+iK 2 ^ 1 (79) 

The latter form of B has been simplified by making the identification x = t\. We see 
that the form of the k dependent factors, or "phase shift functions" appearing in r 2 
enable us to determine the form of the vertex operator relatively easily. This fact has 
been noted elsewhere [11], in connection with the Kadomtsev-Petviashivili equation. 
The form of the operator A is then 

A(k,x) = 1 + e 2 (Er=o K2fe+lt 2fc+i- Ki ) e " 2 (Sfc=o ( 2fc+i) 1 K 2fe+i ei^)^ (go) 

This form of the vertex operator has been cited elsewhere [1,12], and does in fact 
convert an n-bound state r-function into one possessing n + 1 bound states. This is 
most readily seen by direct computation using the identity 

B(K n ,x n )...B( Kl ,x 1 )l=U\\U(^ L ) e 2 fco^ + W^)j (81) 



and comparing, after suitable multiplication by an overall constant and an overall 
exponential factor, with the form of the r-function given in Eq. (|68"1). The n-soliton 
r-function can be generated by repeated action of the vertex operator on ro = 1: 

r n = A(K n ,x n )A(K n „i,x n -i) . ..A(ki,xx)1. (82) 

A particularly compact form of the vertex operator (j80| ) can given if we choose (3 = 
exp(— 2kx) and make use of the fact that B(k,x) 2 = 0. We may write A as 

A(k,P) = e /3B(K ' 0) . (83) 
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In this form, it is apparent that the operator B(k, 0) is a generator of an infinitesimal 
symmetry that maps r-functions onto new r- functions. It has been shown [12] that 
these symmetries are generated by a certain class of infinite-dimensional Lie algebras. 
Another frequently used form of the vertex operator is 

ATM = exp^w) e^(- m+i)K2i+ —J. (84) 



This form of the vertex operator can also be used to add a soliton to a r-f unction: 
we have 

T n = (b+Xfa) + fe-XC-KnJ^Tn-i. (85) 

In this case, it is the ratio b^/b~ that determines the position of the new soliton. The 
equivalence of this construction to that given above can be verified by making use of 
the identity 



x(K n )x(K n - 1 )...x(K 1 )i= n ex pE «4 fc+i w n 



K m ••/,-, V2 



m=l v ' 1=1 



16) 



to construct the n-soliton r-function, and comparing the result, modulo irrelevant 
factors, to that given above. 

In the above discussion we have treated the variables £2fe+i as auxiliary variables 
which are to be set to zero after they have been used to construct a given poten- 
tial. Comparison of Eq. ( flUP with the vertex operator Q8U| ) shows that the auxiliary 
variables t 2 k+i introduced in this section are precisely the KdV time variables em- 
ployed in Sec. IV. Consequently, the above construction can also be used to construct 
multi-soliton solutions of the KdV hierarchy. Defining u(x, t%, t§, . . .) by 

d 2 

u(x,t 3 ,t 5 , ...) = 2— log r(x,t 3 , . . .), (87) 

with t constructed by repeated application of the vertex operator on 1, we see that 
u obeys 

Ut ^ = Tx^ H2k+1 - (88) 

As the time variables run over all possible values, u maps out all of the isospectral 
deformations of the potential V(x) = —u(x)\ t3! t 5 ,...Gxed- 



VI CONCLUSIONS 

We have explored two methods for adding a soliton to a multi-soliton solution of the 
Korteweg-de Vries equation (and related higher-order equations). The method more 
familiar from the standpoint of soliton theory [1] employs a function (the r-function) 
whose logarithm, when differentiated twice, gives the solution (up to a factor). A 
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vertex operator, a function of infinitely many times, acts upon this r-function to add 
a soliton. 

The alternative method for adding a soliton relies upon supersymmetric quantum 
mechanics. In order to add a soliton, one must solve a (first-order) Riccati equation, 
thereby adding one integration constant for each soliton. The physical interpretation 
of this integration constant is the position of the new soliton relative to all the others. 
This position is best visualized at large separations of the individual solitons, which 
occurs at asymptotic values of the times t 2 k+i governing the evolution according to 
the Hamiltonians -f^fc+i- Under such circumstances, the multi-soliton solution closely 
resembles a set of individual one-soliton "lumps," each lump propagating with a speed 
governed by its size. 

We have shown that the two methods for adding solitons are equivalent. Nonethe- 
less, the method of demonstration still seems somewhat roundabout. There are many 
features of single-soliton solutions which suggest that they may be useful in construct- 
ing the vertex operator (|80|), but we have so far been unable to find a more direct 
route to this result. 
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VIII APPENDIX: DERIVATION OF THE TRANS- 
MISSION COEFFICIENT USING SUPER- 
SYMMETRIC QUANTUM MECHANICS 

We wish to show that potentials constructed using the method of Sec. II are reflec- 
tionless and have transmission coefficients T(k) given by Eq. (fL7|) : 

r <*> = nsT£ < 89 > 

These results may be proven by induction on N. 

For, suppose that we have constructed a series of potentials V n , n — 0, 1, 2, . . . , N, 
such that V n has bound states at E = — k\, — k|, . . . , — k\. Using the methods of 
Sec. II, we may construct V n+ \ from V n as follows: We write 

V^ +1) = V n + K 2 n+l = f n+l - (90) 

so that the partner potential V+. n is 

Vi n+1) = V n+1 + K 2 n+1 = f n+1 + f n+ r (91) 
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We therefore find that each potential V n has two equivalent representations: 

Vn — fn+1 ~ fn+1 ~ K n+1 = fn + fn ~~ K n- (92) 

The eigenfunctions of V n can be related to those of V n +i using the operators A and 
A' . Comparing with Sec. II, we find that if ip n is an eigenfunction of V n , then the 
corresponding eigenfunction ip n +i f° r the potential V n+ i is given by 

1p n +l = Al +1 1p n = (— + fn+lj^n- (93) 

We now make the following induction hypotheses: we suppose for some n that (i) 
/„ — > =)=K n as x — > ±00, and that (ii) the transmission coefficient for V n is given by 
Eq. (SS) with N = n. We know from Sec. II that both hypotheses hold for n = 1. It 
follows that V n — > — > ±00, and that / n+ i obeys 

/ n 2 +1 -^ +1 = K + ^ +1 . (94) 

As x — >• ±00, this reduces to f% +1 — f' n+l = K n+1 , which has the solution f n+ i = 
— K n +i tanh K n +i(x — Xo). Consequently / n+ i has the asymptotic forms 

{-n ri+l tanh n n+l (x - x+) -> -« n+ i, x -> +00; 
(95) 
— K n+ i tanh /t n+1 (a; — x„) — > /t n+ i, x — > — 00. 

This proves that our first induction hypothesis holds for all n. To see that the second 
holds, observe that a plane wave solution in the potential V n has the asymptotic form 

e jfcx , x — ► —00; 

^n(x)^<| . (96) 



n?-i ^e ifcx , x -> +00. 



Applying the operator to ^ n gives the corresponding solution for the potential 
Kj.fi. After dividing through by ik + K n+ i, we find that this solution has the form 



e lhx , x — > —00: 



i>n+i{x) -> { (97) 

nS 1 !^^, X-+00, 

which proves the second of our induction hypotheses, as well as Eq. (|B9|). 
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